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Context. Recent observations show that inner discs and rings (IDs and IRs, henceforth) are not preferably found in barred galaxies, a 
fact that points to the relevance of formation mechanisms different to the traditional bar-origin scenario. In contrast, the role of minor 
mergers in the formation of these inner components (ICs), while often invoked, is still poorly understood. 

Aims. We have investigated the capability of minor mergers to trigger the formation of IDs and IRs in spiral galaxies through colli- 
sionless N-body simulations. 

Methods. We have run a battery of minor mergers in which both primary and secondary are modelled as disc-bulge-halo galaxies with 
realistic density ratios. Different orbits and mass ratios have been considered, as well as two different models for the primary galaxy 
(a Sab or Sc). A detailed analysis of the morphology, structure, and kinematics of the ICs resulting from the minor merger has been 
carried out. 

Results. All the simulated minor mergers develop thin ICs out of satellite material, supported by rotation. A wide morphological 
zoo of ICs has been obtained (including IDs, IRs, pseudo-rings, nested IDs, spiral patterns, and combinations of them), but all with 
structural and kinematical properties similar to observations. The sizes of the resulting ICs are comparable to those obtained in real 
cases with the adequate scaling. The existence of the resulting ICs can be deduced through the features that they imprint in the 
isophotal profiles and kinemetric maps of the final remnant, as in many real galaxies. Weak transitory oval distortions appear in the 
remnant center of many cases, but none of them develops a noticeable bar The realistic density ratios used in the present models 
make the satellites to experience more efficient orbital circularization and disruption than in previous studies. Combined with the disc 
resonances induced by the encounter, these processes give place to highly aligned co- and counter-rotating ICs in the remnant centre. 
Conclusions. Minor mergers are an efficient mechanism to form rotationally- supported stellar ICs in spiral galaxies, neither requiring 
strong dissipation nor the development of noticeable bars. The present models indicate that minor mergers can account for the exis- 
tence of pure-stellar old ICs in unbarred galaxies, and suggest that their role must have been crucial in the formation of ICs and much 
more complex than just bar triggering. 

Key words, galaxies: bulges — galaxies: evolution — galaxies: formation — galaxies: interactions — galaxies: kinematics and 
dynamics — galaxies: structure 
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1. Introduction 

Stellar IDs and IRs are found in at least one third of spi- 
ral galaxies (E rwin & Spark e 2002; Falcon-B arroso et a Tl l2006l 
F06 hereafter). Usually masked by the light coming from 
brighter overlapping components, many of them have been 
detected through the peculiar features that these ICs im- 
print in the isophotes or in the ki nematic maps of the 
host galaxy (Ivan den Bosch et alJ 11994: .Kormendv et al.lll994l: 



Lauer et alJ 119951: IScorza et aljTl998l: IScorza & van den Boschl 

1QQSl- "l^r^t^rr^lin T'^ilinrr^r' "OOnrV- "P^cf^fQl' '9001'- "Pr^j/iJ 



1998"; "Koprolin & Zeilinsei^ "200d'; "Rest et al.' '2001'; 'ErwinI 
2004 : Carollo 2004; Liskeretal. 2006; Comeronetal. 2008; 
Morelli et al.ll20iaiMoiseev et al.ll20ia:lSirChenkoll201Q,) . 



Stellar IDs are found in galaxies of all types, although 
they tend to reside in Sa-Sb's. Their radii span from 
some tens of parsecs to ~ 2kpc in diameter. They ex- 
hibit red or blue colors, usually similar to those of the 
host bulges, and harbour stellar popu lations with ages 
ranging from ~ 1 Gyr to ~ lOGyr (iFranx & Tiling worthi 



, 19881: iBender 119901: J van der M arel & Framd 119931 ; 
Cinzano & van de r Marell Il994t ^Prieto et alT l200U 
Carollo et al. 20021 |Erwin & Sparke 2002; Martel et al] 
200 2: .Erwin & Sparke" '2003; Kormendv & Kennicutl .2004j: 



Mo relli etal.1 12 004: Athanassoula 2005; Cappellari et al. 1 120071: 
lEmsellem et al.l l2007i: iPeletier et al.i i2007l) . Inner star-forming 
IRs and pseudo-ri ngs are detected in at least one fifth of 



all disc galaxies ([Knapenl l2005l : iMoiseev & BizvaevI l2009l : 
Ide Lapparent et al .11201 1 "). However, accounting for the episodic 
nature of their star formation histories and their long-lived 
stable configurations, the total fraction of IRs including pure 
s tellar (non star-forming) ones is expe cted to be much higher 
(iMazzuca et al.ll2006l:[Sarzi et al.ll2007l) . 

These stellar ICs could be primordial features in the 
galaxy centres, established at an early epo ch of rapid , 
cold accretion flows at high redshift (Bedregal et al.l |2006|). 
However, the dramatic decline of the star formation surface 
density in the discs during the last ~ 8 Gyr and the correla- 
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tions found between the nuclear and disc global properties 
of spiral galaxies evidence an intertwined star formation 
history for bulges, discs, and I Cs that spreads through 
time (Peletier & Balcells 1996: IPominguez-Palmero et alJ 
I2QQ8I: iDominguez-Palmero & Balcellsl l2QQ8l) . Attending 
to the high detection fraction of stellar IDs and IRs in 
barred galaxies, their origin has be en traditiona ll y as- 

sociated to bar patterns in discs (jCombes et alj Il992t 

Arnaboldi et al. 1995a 



'95i; Ivan den Bosch & Emsellem 1998; 
m iRegar 



Erwin & Sparke 20021: iRegan & Teubea .2003: .Bokeretal.i 



20Q8[ Ide Lorenz o-Caceres et al.l 12008^ IComeron et al.l l2010h . 
This has been corroborated by numerical simulations, that 
have shown that the dynamical resonances induced by bars 
give place to the rings surrounding them easily. Moreover, 
the bars induce strong gas inflows to the galaxy centre, which 
settle in cir cular orbits, in ducing star formation that can resu lt 
in an ID ( Norman etaP 119961: lAthanassoula et alJ l2009bllal: 
IComer6n"et al. 2010). Other bar-related scenarios propose that 
IDs can also be the relics of diluted nested bars or of episodes 
of starbursts in rings that have shrunk in radius to the galaxy 
centre over ti me (Combes et al. 1992; Regan & Teuben 2003; 
ISirChenko & Smirnoval l2010). 

Nevertheless, recent observations indicate that stellar IDs 
and IRs are not preferably harboured by barred galaxies, being 
at least as frequent in non-barred early-type hosts as in barred 
ones dE msellem et al . 2004; Fal con-Ba rroso et al. 2004; Knapen 
| 2005blKnapen et al.l l2006l : ISarzi et al] l2006l: fcomeron et al.l 
12010'). IRs are easily developed by unsta ble gas-rich discs in nu- 
merical simulations ( Aumer et al .1 120 1 Oh . b ut IDs are structure s 
more difficult to reproduce spontaneously (' Thakar et al.l 1 19971) . 
So, other mechanisms capable of triggering disc resonances be- 
sides bars are required to explain the existence of stellar IDs and 
IRs in unbarred galaxies. One of the main alternatives is gas 
infall, but simulations indica te that this mechanisn i is efficient 
producing IRs, but not IDs (IThakar & Rvdenlll998l: iNaab et al.1 
120061) . 

Another possible driver is merging, as suggested by the 
existence of many ICs exhibiting strong misalignments or 
distorted morphologies and/or kinematics with respect to the 
host disc (see, e.g.. lOkumura et al.l 11994 I A rnabold i et al.l 



1995bl: iButa & C ombes' '1996': 'Barnes & Hemauist 



Knapen etal.1 120 04: Reshetnikov et al. 2005; Mazzuca et al. 



2006; Sil'chenko & Moiseev 2006; Chilingarian et al. 2009; 
f aundez-Abans et al. 2009; Moiseev & Bizvaev 2009; 
iBrosch et al.l l2010h . This scenario is supported by numeri- 
cal simulations, which have proven that major mergers can 
drive the formation of kinematicall y-decoupled ICs analogous 
to the o nes found in E-SO galaxies ("Hemauist & Barnes" 1991'; 



Balceffs & Gonzalez 1998; Bendo & Barnes 2000; Barnes 2001; 



Jesseit et al.i 2007: Di Mat teo et al.ll2008l) . Considering that mas- 



sive E-SO galaxies have experienced at le ast one major merger 
since z_~ 1 ( Gonzalez-Garcia et al. 20091; iLopez-Sanjuan et all 



120091: Imiche -Moral et al. 2010a,b), this means that the existence 
of ICs in E-SO's can be satisfactorily explained through this 
process. However, they can not account for the ICs found in 
unbarred Sa-Sb galaxies, as the remnants resulting from major 
merge rs are E-SO's for typical gas amounts (Barnes & Hemauist 
Il996l: iBoumaud et al.ll20Q4 l2005bl: iNaab & Burkertll2003>) . ' 

The straightforward altemative to major mergers is minor 
merging. Minor mergers induce dynamical resonances and oval 
distortions in discs easily and imprint a smooth growth to the 
pre-ex isting galax y bulge (see Bertola et al. 1999; Pizzell a et al.l 
20021: lE Hche-Moral et al.1 120061: ISil' chenko & Moiseevll2006 



Brosch et al...201Ql) . Although the eff'ects of minor mergers on 



bulge and disc gro wth and on the satel l ites have been exten- 
sively studied (e.g., Hopkins et al. 2010; Henriaues & Thomad 
2010; Tapia et al. 2010a b; Bartoskova et al. 2011; Ebrova et al.] 
I2OII ), little attention has been devoted to study specifically their 
ability to induce the formation of dynamically-cold ICs in galax- 
ies. 

One of the earlie r studies dealing with this topic was 
carried out by Elme green et al.l (Il992), who performed nu- 
merical simulations of ring-companion interactions to anal- 
yse the effects of the minor merger onto pre-existing outer 
rings. Later, Thakar and collaborators tested the formation 
of counter-rotating discs in spiral galaxi es through retrograde 
mergers of gas-rich dwarfs onto discs (Thakar & Rydenlll996[ 
iThaikar et al.lll997l:lThakar & RvdenllT998h . Thev found that, al- 
though counter-rotating thin gaseous discs were formed during 
the minor merger, the obtained sizes were only comparable to the 
biggest observational cases. The formed IDs did not have expo- 
nential profiles either and were highly unstable, quickly deriving 
into IRs. This led these authors to conclude that, in order to form 
a normal counter-rotating disc, "there must be either little or no 
pre-existing prograde gas in the primary galaxy, or its dissipative 
inffuence must be off'set by signific ant star formation activity". 
More recently. 'A guerri etaP (120011) performed collisionless N- 
body simulations to test the growth of bulges after the accretion 
of dense spheroidal satellites. Their undisrupted dense satellite 
cores sank to the galaxy centre, giving place to kinematically- 
decoupled components in the remnants, but not supported by ro- 
tation^ 

lEhche-Moral et al.l (I2006L EM06 hereafter) studied the ef- 
fects of the accretion onto disc galaxies using satellites which 
themselves comprised a disc, a bulge, and a dark halo. More im- 
portantly, the relative densities of primary and secondary were 
made to be realistic by imposing that both models lie on the 
Tully-Fisher relation. This was an improvement over previous 
studies of mass buildup via accretion, given that the disruption of 
the satellite, the radius of deposition of satellite material, and the 
dynamical heating of the primary, all depend on the tidal fields, 
which scale with the relative densities of the two galaxies. EM06 
reported the formation of dynamically-cold stellar structures in 
the centre of remnants, made out of dismpted satellite material. 
However, that paper was centred on the bulge growth driven by 
the minor merger, and hence no analysis or description of the 
ICs was performed. We have therefore extended the simulations 
of EM06, sampling different orbits and initial conditions and us- 
ing ~ 3 times more particles, to carry out an exhaustive study 
of the morphology and kinematics of the ICs resulting from dif- 
ferent minor mergers. We have simulated collisionless cases, as 
recent studies have demonstrated that dissipative components 
are not decisive in the formation and shaping of kinematically- 
decoupled ICs through major mergers (although gas makes them 
more axisymmetric and can induce recent star formation in them, 
see Jesseit et al. 2007; Di Matteo et"aDl2008h . 

The present models demonstrate the capability of minor 
mergers to induce the formation of a wide morphological zoo 
of thin stellar ICs in the galaxy centres, with stmctural and kine- 
matical properties analogous to those harboured by real spiral 
galaxies. The novelty of these models is two-fold: first, all the re- 
sulting ICs are made out of dismpted satellite material (whereas 
they come from resonances in the parent disc in previous studies) 
and, secondly, dissipative effects and strong bars are not essen- 
tial to form these dynamically-cold ICs. 

The paper is structured as follows. We brieffy describe the 
models in ^ The formed ICs are analysed geometrically, photo- 
metrically, and kinematically in ^ In ^ we perform a qualita- 
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Table 1. Number of particles used in the models 







Number of 


Particles (/ 10^) 




Experiment type 


Total 


Dl 


Bl 


HI 


D2 


B2 


H2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Big bulge 


185 


40 


10 


90 


10 


5 


30 


Small bulge 


415 


60 


10 


300 









Columns: (!) Experiment type depending on the used primary galaxy 
model (big or small bulge). (2) Total particle number. (3) Number 
of primary disc particles. (4) Number of primary bulge particles. (5) 
Number of primary halo particles. (6) Number of satellite disc parti- 
cles. (7) Number of satellite bulge particles. (8) Number of satellite 
halo particles. 

tive comparison of the ICs resulting in the models with those de- 
tected in real spiral galaxies. Model limitations are commented 
in ^ The discussion can be found in ^ A brief summary of 
the results and some conclusions are finally addressed in ^ 
The physical magnitudes are provided in units of the simula- 
tion throughout the paper, although scaling to real galaxies is 
straightforward assuming the scalings described in ^ 

2. The models 

We have extended the set of collisionless N-body simulations 
of minor mergers onto disc galaxies described in EM06, now 
using longer pericenters and diff'erent initial disc galaxies. The 
outcome of the satellite material of a final set of 12 collisionless 
models has been analysed (six of them come from EM06). Ten of 
these experiments were run using a disc galaxy with a prominent 
bulge as primary galaxy (equivalent to a Sa-Sb galaxy), while 
in another two a primary with a smaller bulge (similar to a Sc) 
was used to investigate the influence of the primary bulge in the 
outcome of the accretion. 

All the galaxies in the simulations (primary galaxies and 
satellites) have an initial bulge-disc-halo structure. The pri- 
mary galaxy models were built using the GalactlCS code 
(^uijken & Dubinski 1995), including an exponential disc com- 
ponent ( Shu[1969|), a King bulge ( Kins 196 ^, and a dark hal o 
built following an Evans profile ( Kuijken & Dubinskil Il994l) . 
The discs of both primary galaxy models follow an exponen- 
tial surface density profile both radially and vertically, and were 
allowed to relax in isolation for about 10 disc dynamical times 
prior to placing them in orbit for the merger simulations. No rel- 
evant resonant structures appear in the disc. The primary galaxy 
with a big bulge matches the Milky Way (MW) when the units 
of length, velocity, and mass are R = 4.5 kpc, v = 220 km s"^ 
and M = 5. 1 X lO^^M©, respectively. In this case, the correspond- 
ing time unit is 20.5 Myr. The primary galaxy with a small bulge 
matches NGC 253 using the following units of length, velocity, 
and mass: R = 6.8 kpc, v = 510km s"\ andM = 2.6 x lO^^M©, 
implying a time unit of 1 1.7 Myr. These values, especially when 
using an appropriate M\^^IL ratio, yield mass-to-light ratio val- 
ues close to observations (M/L -10). Tables[T]and[2l summarize 
the main characteristics of the two relaxed galaxy models used 
to represent the primary galaxy. For the remainder of the article 
we use model units in all the figures, a conversion to physical 
units is easily achieved through the above correspondences. 

Satellites are scaled replicas of the primary galaxy model 
with a big bulge in all the experiments. A physically-motivated 
size-mass scaling was used to ensure that the primary-to- satellite 
density ratios are rea listic, forcing both g alaxies to obey the 
Tully-Fisher relation (iTullv & Fished 1 19771 consult EM06). In 



Table[3l we list the characteristic masses and sizes of the satel- 
lites for each merger experiment. The luminous mass ratios be- 
tween the primary galaxy and the satellite considered have been 
1:6,1:9, and 1 : 1 8. In the models with big primary bulges, the pri- 
mary and satellite are scaled replicas, so the luminous mass ratio 
of the encounter is equivalent to its total mass ratio. However, 
as commented above, the satellites in the models with small pri- 
mary bulges are scaled versions of the big primary bulge model, 
in order to analyse the influence of the primary bulge structure 
in the outcome of the accretion. Satellites were also relaxed in 
isolation prior to the merging simulation. 

The initial separation of both galaxies was 15 primary disc 
scale-lengths in all the experiments. In order to avoid a perfect 
spin-orbit coupling, the initial inclinations between the orbital 
plane and the galactic planes of the primary and the satellite 
galaxies were fixed to 30° in direct orbits and to 150° in ret- 
rograde orbits, respectively. All satellites have also an azimuthal 
angle = 90°. For each satellite mass, pericenter distance, and 
model of the primary galaxy considered, a direct and a retro- 
grade orbit have been computed. Initial orbits were elliptical 
with pericenters equal to one or eight disc scale lengths, de- 
pending on whether we were running a short or a long pericen- 
ter orbit. Relative velocities at the first pericenter passage oscil- 
late between 440-650 km/s, considering the scalings commented 
above. 

In Table[3l we include the orbital parameters of each merg- 
ing experiment. As satellite models are scaled down versions of 
the primary model with a big bulge and exhibit diff'erent mass 
ratios having the same number of particles, this renders dif- 
ferent mass particles for the diff'erent components in each ex- 
periment. The largest mass contrast, and thus the higher two- 
body errors, are expected in the models with a total mass ratio 
of 1:18, where a primary halo particle is 10 times more mas- 
sive than the bulge particles of the satellite. Such extreme ra- 
tio is well below the limits explored earlier in diff'erent sim- 
ulations using a similar set of initial models, so the kinemat- 
ics of the inner regions of the remnant a re much better sam- 
pled in the present models than previously ('Ba lcells & G onzalei 
1998; Gonzalez-Garcia & BalceUs 2005, EM06). The diflTerence 
in the number of particles for the haloes hosting a big and 
a small bulge is two folded. On one hand, a bulge-less disc 
mod el stable to bar distortions requires a concentrated halo (see, 
e.g., Gonzalez-Garcia & BalceUs 2005). On the other hand, the 
masses of the halo and bulge particles need to be of the same 
order for ensuring low two-body errors. These two constraints 
required the models with small bulge to have such a high num- 
ber of halo particles. As we are interested in the physics of the 
inner regions of the remnants, we intended to balance between 
the accuracy gained by using a higher number of particles and 
the economy of resources in building new stable models for each 
satellite. 

The evolution of the new models was compu ted using 
GADGET-2 code ( Springel et al .ll2001l:ISpringdll2Q05h . We used 
a softening of s = 0.02 in model units and an opening angle 
of ^ = 0.6. Considering this tolerance parameter and apply- 
ing quadrupole-moment corrections, the code computes forces 
within 1 % of those given by a direct summation and preserves 
total energy better than 0.1%. We evolved all models for ~4 halo 
crossing times beyond full merger to allow the final remnants to 
reach a quasi-equilibrium state with a good conservation of en- 
ergy in all runs. Times to full merger and total run times of each 
experiment are also indicated in Table[3]in simulation units. The 
initial B/D ratios of the main galaxy increase after the merger to 
B/D ^ 0.6 in the cases of big bulges and B/D ^ 0.2 for small 
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Table 2. Initial parameters of the primary galaxy models 



Experiment type 




Mb,i/Mb,i 


A^Dark,l/>tL,l 


hD,i 


rB,l/Vl 


^95%,l/^D,l 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Big primary bulge 


6.40 


0.51 


4.16 


1.00 


0.20 


3.7 


0.11 


Small primary bulge 


0.99 


0.08 


6.98 


0.39 


0.16 


3.8 


0.06 



Columns: (!) Experiment type depending on the primary galaxy model used (big or small bulge). (2) Total primary mass (simulation units). (3) 
Bulge-to-disc mass ratio of primary galaxy. (4) Dark- to-luminous mass ratio of primary galaxy. (5) Radial disc scale-length of the primary galaxy. 
(6) Ratio between the effective radius of primary bulge and primary disc scale-length. (7) Ratio between the radius of the shell containing 95% 
of luminous material in the primary galaxy and its disc scale-length. (8) Ratio between the vertical and radial disc scale-lengths of the primary 
galaxy. 



bulges. This means that the remnants of the experiments start- 
ing with a Sb primary galaxy have become SO-Sa's, while the Sc 
primary galaxies have been transformed into Sb-Sbc's after the 
minor merger (see Graham 2001). 

We will refer to each model throughout the paper according 
to the following code: MmP[l/s][D/R][b/s], where m indicates 
the bulge-to- satellite mass ratio (m = 6, 9, or 18 for models with 
luminous mass ratios equal to 1:6, 1:9, 1:18, respectively), "PI" 
indicates long pericenter and "Ps" short pericenter, "D" or "R" 
describes the orbit ("D" for direct and "R" for retrograde), and 
the final "b" or "s" letter indicates if the primary galaxy had a 
big or a small bulge (see TableO. 



3. Results 

In EM06, we reported the formation of inner dynamically-cold 
components in minor merger experiments, but we did not anal- 
yse the resulting ICs neither structurally nor kinematically. In the 
next sections, a detailed analysis of the co- and counter-rotating 
ICs formed in the models described in ^is performed. 

3.1. Formation of ICs 

In Fig. [T] we show the disruption experienced by the satellite in 
model M18PsDb (model i in Table[3l) as an example of the time 
evolution of the luminous surface density of the satellite ma- 
terial during the last moments of the encounter. Although the 
diff'erent initial conditions aff'ect to the global shape, size, and 
even the number of components in these inner structures, the fi- 
nal structure in all the experiments resembles a central vertically 
thin torus or disc (depending on whether the satellite material 
reaches the remnant centre or not), embeded into a more ex- 
tended flat component, similar to a more extended disc. In the 
model plotted in the figure, the toroidal structure corresponds to 
the central ring visible at 7? ~ 0.5, while the outer disc corre- 
sponds to the low density structure that extends up to 7? ~ 1-1.5. 
The morphologies of the ICs formed in each experiment are de- 
scribed in detail in ^3.2.31 

Figure[T] also shows that the outer and inner regions of the 
structure resulting from the satellite disruption are built up at 
diff'erent epochs during the minor merger. The less bound parti- 
cles of the satellite (i.e., those from the disc) are disrupted earlier 
in the interaction. During the first pericenter cross of the orbit, its 
outer shells are removed from the satellite by the primary galaxy 
tidal field, giving place to the outer structure of the formed IC 
(as observed in the first half of snapshots in Fig.[T|). The core of 
the satellite takes more time to experience a noticeable disrup- 
tion. Its material is deposited at inner radii during the last stages 
of the encounter. 



3.2. Geometrical and kinematical characterization of ICs 

In this section, we describe the geometry and structure of these 
ICs through their radial and vertical surface brightness profiles, 
and analyse their misalignment with respect to the galaxy plane 
of the final remnant and their kinematics. 

3.2.1 . Radial distribution of ICs in the remnants 

FigureO shows the radial surface density profiles of luminous 
material initially belonging to the satellite, to its disc, and to its 
bulge in the final remnants of all the experiments. In general, the 
inner structure is basically composed by satellite bulge particles 
(typically, at R < /zd,i(^ = 0)), while the outer parts are con- 
trolled by the surface density of satellite disc particles. However, 
there is no sharp transition between the two radial ranges where 
each satellite component dominates. 

We have fitted one exponential disc or a combination of two 
nested ones to the surface density profiles of the satellite stellar 
material in each remnant (represented by two solid black straight 
lines in each panel). We assume that the satellite material de- 
scribed by each exponential radial profile defines an independent 
structure inside the global IC formed in the remnant center out of 
disrupted satellite material. In the fits, we have rejected the cen- 
tral regions of the density profiles to avoid the particular features 
exhibited by each distribution in the centre, such as undisrupted 
satellite cores (as in M6Ps[D/R]s, models e and f), holes (as in 
M18P[l/s]Db, models i and k), and smooth-length eff'ects. The 
resulting fits and their residuals are plotted in the corresponding 
panels and sub-panels of Fig.[2l 

Most of the resulting ICs can be well- approximated by 
one exponential radial profile or by the addition of two ones. 
However, notice that the independent ICs identified by each fit- 
ted exponential radial profile are shaped with material of both 
satellite components: from the bulge and the disc. Only in the 
experiments with a small primary bulge, the satellite bulge ends 
undisrupted in the remnant centre (see the profiles corresponding 
to the satellite bulge particles in M6PsDs and M6PsRs, frames e 
and f in the figure). In these models, the outer shells of satellite 
material configure an extended ID structure hosting the undis- 
rupted satellite core. 

In order to delimit the radial extent of the ICs characterized 
by an unique exponential profile, we have assumed that, in the 
case that only one exponential profile is required to explain the 
whole radial structure, this IC extends up to the radius at which 
the fitted surface density is equal to 1/10 of its central value as 
extrapolated from the fit. In the case that the structure is bet- 
ter described by two nested exponential profiles, the innermost 
one is considered to extend up to the radius where the outer one 
starts to dominate the global fit. The outer one will extend from 
this radius up to that at which its surface density drops to 1/10 
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Table 3. Orbital and scaling parameters of each merger experiment 



Model Code 


Code in EM06 


X XXlllClX y ±J KA-L^K^ 




-^Sat/'^Prim 


-'^pencenter/ '^D,l 


'^Prim 


^fuU merger 


* total 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(a) M6 Ps Db 


M2TF35D 


Bis (h) 


1-6 rM6^ 


0.46 


73 (VO 


30 (D) 


-72 


100 


(b) M6 Ps Rb 


M2R 


Big (b) 


1:6 (M6) 


0.46 


0.73 (Ps) 


150 (R) 


-80 


100 


(c) M6 PI Db 




Big (b) 


1:6 (M6) 


0.46 


8.25 (PI) 


30 (D) 


-93 


144 


(d) M6 PI Rb 




Big (b) 


1:6 (M6) 


0.46 


8.25 (PI) 


150 (R) 


-110 


144 


(e) M6 Ps Ds 





Small (s) 


1:6 (M6) 


0.25 


0.87 (Ps) 


30 (D) 


-40 


62 


(f) M6 Ps Rs 




Small (s) 


1:6 (M6) 


0.25 


0.87 (Ps) 


150 (R) 


-44 


72 


(g) M9 Ps Db 


M3TF35D 


Big (b) 


1:9 (M9) 


0.39 


0.79 (Ps) 


30 (D) 


-80 


100 


(h) M9 Ps Rb 


M3R 


Big (b) 


1:9 (M9) 


0.39 


0.79 (Ps) 


150 (R) 


-87 


100 


(i)M18Ps Db 


M6TF35D 


Big (b) 


1:18 (M18) 


0.28 


0.86 (Ps) 


30 (D) 


-116 


122 


G)M18Ps Rb 


M6R 


Big (b) 


1:18 (M18) 


0.28 


0.86 (Ps) 


150 (R) 


-142 


154 


(k)M18Pl Db 




Big (b) 


1:18 (M18) 


0.28 


8.19 (PI) 


30 (D) 


-225 


260 


(1)M18P1 Rb 




Big (b) 


1:18 (M18) 


0.28 


8.19 (PI) 


150 (R) 


-285 


340 



Columns'. (1) Model code: MmP[l/s][D/R][b/s], where m indicates the bulge-to- satellite mass ratio (m = 6, 9, or 18 for models with luminous mass 
ratios equal to 1:6, 1:9, 1:18, respectively), "PI" refers to long pericenter and "Ps" to short pericenter, "D" indicates direct orbits and "R" retrograde 
orbit), and the final letter ("b" or "s") indicates if the primary galaxy had a big or a small bulge. The letter in parentheses helps to identify each 
model quickly in the forthcoming figures. (2) Model code in EM06, for those models that were already presented in that paper. (3) Primary galaxy 
model used in the experiment (big or small primary bulge, see Table[T]). (4) Luminous mass ratio between satellite and primary galaxy. (5) Ratio 
between the luminous half-mass radii of the satellite and the primary galaxy. (6) First pericenter distance of the orbit, in units of the primary disc 
scale-length. (7) Initial angle between the orbital momentum and the primary disc spin. This angle determines if the orbit is prograde (direct) or 
retrograde. (8) Approximate time of full merger, in simulation units. (9) Total run time of each experiment, in simulation units. 



of its central value. The radial extent of each radial component 
resulting from these criteria are marked in Fig.E] 

Summarizing, all the simulated minor mergers give place to 
complex extended ICs in the remnants made out of disrupted 
satellite material, with radial surface density profiles that can be 
well-described by one or two nested exponential profiles. 

3.2.2. Vertical distribution of ICs in the remnants 

An IC with a radial exponential profile can correspond to an ID 
or to a bulge with Sersic index n - \ (i.e., a pseudo-bulge). 
The diflTerence between all these components arise in the vertical 
thickness of the component, as compared to its o wn radial scale- 
length and to the one of the host disc. Following Sil'chenkoet 
(|2Q1 1), a ratio of scale-length to scale-height of about 3 is a rea- 
sonable frontier between spheroids and discs. So, we have esti- 
mated the ratio between the vertical and the radial scale-lengths 
of the formed ICs in each remnant, as well as the ratio of their 
scale-heights to the radial scale-length of the remnant discs, in 
order to identify them as IDs (thin or thick) or pseudo-bulges. 

Figure[3] presents the characteristic vertical scale-lengths /Zz 
at diflTerent radial positions of the disc, bulge, and all luminous 
material originally belonging to the satellite in the final rem- 
nants. The scale-lengths have been derived from exponential fits 
to the vertical density profiles at each radius. Only the radial po- 
sitions with enough particles to ensure a vertical density profile 
of S/N > 50 have been considered (typically, R < 2/zD,primary). 
The fitting errors are on average < 10% for each estimate. The 
final scale-length of the disc remnant is marked in each panel as 
a reference (horizontal dotted lines). We have also indicated the 
radial extent of each IC, according to the definition adopted in 

Notice that the IC made out of satellite material presents 
scale-heights smaller than the radial scale-length of the original 
primary disc along the whole considered radial range in all the 
models, although it is vertically wider than the structure made 
out of primary disc material at all radii in all the remnants (com- 
pare asterisks and squares in Fig.|3j. So, although the primary 



disc has been heated by the satellite accretion, it is thin as com- 
pared to the IC. Therefore, the satellite material contributes to a 
structure thicker than the main remnant disc. This result is co- 
herent with the structures of disc galaxies resulting from ACDM 
cosmological simulations, where the disrupted sa tellites con- 
tribute to the buildup of the thick disc of the galaxy (lAbadi et al.l 
[2003). In general, the scale-height of the structure made of dis- 
rupted satellite material increases with radial position, meaning 
that these structures are flared (see asterisks in the figure) . 

The material of the disrupted satellite bulge tends to exhibit 
lower scale-heights than that of the disrupted satellite disc, but 
the diflTerent disruption epochs of both components seems to not 
aflTect to their final vertical distributions at a given radial posi- 
tion, as the scale-heights of both distributions are similar at each 
radius (compare red diamonds and blue triangles in the figure). 

In Fig.O we also show the ratio between the vertical and 
radial scale-lengths of each IC (numbers in black characters in 
each frame). Nearly half the innermost ICs formed in the models 
with big primary bulges are thin, as they exhibit ratios typically 
below ~ 0.4. Thi s means that these inne rmost components are 
IDs, according to lSil'chenko et al criterion. Other cases 

exceed the limiting value proposed by these authors (as models 
M6P[s/l]Rb, panels b and d in the figure), meaning that these 
ICs would be classified as pseudo-bulges in the case that the 
original primary galaxies lacked of a large central component. 
As this is not the case (they had a massive central bulge), the 
ICs formed in these models are embeded into the pre-existing 
bulge or the galaxy thick disc. We also have some question- 
able cases with vertical-to-radial scale-lengths near the limiting 
value for distinguishing between IDs and pseudobulges, as mod- 
els M[9/18]PsRb (panels h and j in the figure). 

The ratios of the scale-heights of these innermost compo- 
nents to the radial scale-length of the remnant disc indicate that 
the majority of these innermost ICs exhibit ratios below 0.4 (see 
the numbers in blue characters at the top left of each panel in the 
figure). This means that these structures are buried in the final 
remnant disc. The models with small primary bulges (panels e 
and f in the figure) exhibit high scale-heights in the center, be- 
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Fig. 1. Time evolution of the luminous surface density of the satellite material during the last moments of experiment MlSPsDb 
(experiment "i" in TableO. Face-on and edge-on views centred in the initial primary galaxy are plotted (left and right panels, 
respectively). Snapshots corresponding to times from 101 to 120 are shown from up-to-down and left- to-right of each figure, using 
a time step equal to one. If the primary galaxy is scaled to the MW (see ©, the total time period represented corresponds to 
~ 0.5 Gyr. A rainbow colour palette is used to represent diff'erent surface density levels in logarithmic scale, with redder colors 
indicating higher values. Spatial scales in both axes are provided in simulation units. The disrupted satellite material is finally 
deposited in the remnant centre forming an IR (in yellow) embeded into a more extended ID (in green). The IR rotates in the same 
direction as the primary disc material. 



cause the satellite core sinks to the remnant center without expe- 
riencing disruption. None of these two models develops a pseu- 
dobulge out of disrupted satellite material, because the undis- 
rupted satellite core produces a central peak in the radial surface 
density profile with a Sersic index n - 2 (see the corresponding 
panels in Fig.[2l). Some of the innermost ICs also present a cen- 
tral hole in their radial surface density profiles (see panels d, i, 
and 1 in Fig.[2l). So, considering their low thickness, they must be 
IRs instead of IDs (this is confirmed by surface density maps in 
E231). 

The outer ICs are thick in general, following the 
[Sir chenko et"an criterion, so we can conclude that these struc- 
tures basically shape the thick disc in the remnants. 

Therefore, the minor mergers give place to ICs in the rem- 
nant centers out of satellite material that can be described as sim- 
ple IDs, nested IDs, IRs hosted by IDs, and bulges hosted by IDs. 
Although their thickness depends on the case, the central regions 
of the ICs are generally embeded in the remnant disc and are rel- 
atively thin, while the outer ones constitute a thick flared discy 
structure surrounding the main remnant disc. The global struc- 
tural characteristics of each IC as derived from their radial and 
vertical density profiles in Figs.|2K3]are summarized in Tabled 



3.2.3. Morphology of ICs 

Figures|4l[5]show the morphology of the ICs formed in each rem- 
nant. In order to provide a 3D description of the ICs structure, 
we plot the face-on and edge-on surface density maps of the 
disrupted stellar satellite material in all the remnants. The fig- 
ure shows that the satellite disruption in our experiments has 
given place to a wide zoo of ICs, all made out of disrupted satel- 
lite material, and confirms the global morphology derived from 
the radial and vertical surface density profiles of Figs. |2]|3] (see 
Tabled. 

Moreover, the maps reveal the existence of substructures 
than make the morphology of the IC even more complex, such 
as ring relics (MlSPlDb, panel k in Fig.O, pseudo-rings (i.e., a 
non-closed ring-like structure like the one developed in exper- 
iment MlSPlRb, panel 1 in Fig.O, and spiral arms (MlSPsDb, 
panel i in Fig.O. The twin clumps observed in some ICs could 
be tracing weak bars or relics of oval distortions (see panels a, 
c, and d in Fig.|4]and panels h-i in Fig.O. The global ellipticity 
and PA profiles of the final remnants seem to support the bar- 
related origin of some of these structures (see ^3.3.11) . However, 
none of the remnants has developed a clear nuclear bar in the 
IC, despite of having IRs and pseudo-rings. We summarize the 
substructures found in each case in Col. 3 of Table|4] 
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Our models demonstrate that minor mergers onto disc galax- 
ies can give place to dynamically-cold thin ICs, with varied mor- 
phologies ranging from nuclear bars embeded in nested IDs to 
pseudo-rings, covering all different kinds of inner structures. 

3.2.4. Bars and oval distortions induced by the minor merger 

Numerical simulations have shown that tidal interactions easily 
produce bars and non-axisymmetric distortions (such as ovals) 
in thin discs through gravitational torques (see jB ournaud et alj 
2QQ5al: >Mastropietro et al] 120051: lAguerri & Gonz alez-Garcial 
20091) . The time evolution of the density maps of the primary 



disc material reveals that our simulations are not exceptions: 
elongated distortions in privileged directions, defined by the or- 
bit of impact, are induced in the primary disc by the minor 
merger (see Fig.O. Some of these transient bar-like distortions 
remain until the end of the simulation, but extremely weakened 
(as in the case of the figure), but most of them dissolve during 
the last stages of the remnant relaxation. 

The strongest and longest-lived central oval distortions in- 
duced by the minor merger in the primary disc appear in the 
cases with a small primary bulge, corroborating many previ- 
ous studies that have posed that a high mass concentration in 
the galaxy centre tends to stabilize the disc, prev enting self 
gravit y and thus bar formati on (see , e.g.. iPfenniger & Normad 
19901: iB oumaud & Combes' "2002"; 'A thanassoulaet al.i ,2001 
EHche- Moral et al. 2006; Cox et al. 20081). AH the transitory bar- 
like or oval patterns formed in the primary discs rotate in the 
same direction as the primary disc stars even in the retrograde 
mergers. 

The final distribution of primary disc material does not show 
distinct ICs, except in experiment M9PsRb (see Fig.[7]), where 
the primary disc material develops a pseudo-ring at the location 
of the twin clumps observed in the IC made out of satellite mate- 
rial (see panel h in Fig.O. Therefore, we are led to conclude that 
the minor mergers do not drive the formation of IDs or IRs made 
out of primary material in our experiments, but force the intro- 
duction of primary disc material to the centre instead, making 
the final galaxy bulge to grow larger (as reported in EM06). 

The resonances that the minor merger induce in the primary 
disc couple with the satellite disruption, triggering the forma- 
tion of ICs, such as IDs and IRs (see panels d in Fig.[4l and i 
in Fig.O. Nevertheless, no clear nuclear bars are developed in 
the IC formed out of disrupted satellite material either. So, our 
models prove that minor mergers can give place to the formation 
of IRs without requiring the development of noticeable bars. 



3.2.5. Alignment of ICs in the remnants 

We have analysed the alignment of the diff'erent inner structures 
formed in the merger with respect to the global remnant struc- 
ture. This analysis is restricted by the fact that our models ex- 
plore only two orbit inclinations (30° and 150°), hence little can 
be said on the eff'ect of diff'erent inclinations in the alignment of 
the resulting ICs. In Fig. [HI we plot the angles between the an- 
gular momenta of the stars originally belonging to the satellite 
bulge and disc in the final remnant and its total angular momen- 
tum (diamonds and triangles, respectively). This figure shows 
that there is a nearly co-planarity between the IC formed out of 
satellite bulge material and the final galactic plane of the rem- 
nant in all the models (notice that their angles are ~ 0° in the 
direct orbits and ~ 180° in the retrograde ones). The alignment 
of the structure formed out of satellite disc material is slightly 





2 f b) M6PsRb: <hJ/h,(ICs)=0.97,0.^01 0.40,1.04 - 
. . ^ . 1 

q ■ - I 




a) M6PsDb: <hJ/hr(iCs)=o.38,o.82 0.05,0.44 




<hjj/hr(ICs)=0.35,0.26 0.03,0.09 " 



d) |vi6PlRb: <hj/i|,(ics)= 1.93,1. 05 0.16,0.40 
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Fig. 3. Vertical exponential scale-lengths of material in the final 
remnants as a function of radius. The plotted physical magni- 
tudes are provided in simulation units. The letters identifying 
each panel are those used in Table[3]for each model. Triangles: 
Material originally from the satellite disc. Diamonds: From the 
satellite bulge. Asterisks: Stellar material originally in the satel- 
lite. Squares: From the original primary disc. Horizontal dotted 
line: Scale-length of the final remnant disc, marked as a refer- 
ence. Vertical dashed lines: Radial extent of each IC character- 
ized by an unique exponential profile, according to Fig.|2l The 
ratio of the vertical to the radial scale-lengths of each IC is in- 
dicated in black characters (one number per radial exponential 
ICs). The ratio of the scale-height of each IC to the radial scale- 
length of the final remnant disc is also indicated in blue charac- 
ters at the top left of each panel. 



lower. The near perfect alignment indicates that the plane of the 
orbit of the satellite has precessed to that of the primary disc 
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Table 4. Characteristics of the ICs formed in the remnants 



Model 


Global structure 


Substructure 


Rotation 


Detectable features 


(1) 






(4) 




(a) M6 Ps Db 


Nested IDs (thin and thick) 


Clumps 


Co-rot. 




\0) IVID r S ivD 


iNesieQ ius i^Dom inicK^ 


Oval distortion 


Counter- rot. 


/z4 > at the ID location. 


(c) M6 PI Db 


Np<^tprl IDs rbnth thin") 




Co-rot. 


^ anH PA trpnrl<s of wpak oval rli<stnrtinn 


(d) M6 PI Rb 


IR+ID (thick) 


Clumps 


Counter-rot. 


6 and PA trends of weak oval distortion. 










Ring sections in and h4 maps. 










Dumbbell structure in <t map. 


(e) M6 Ps Ds 


Bulge+ID (thick) 


— 


Co-rot. 


— 


(f) M6 Ps Rs 


Bulge+ID (thick) 


— 


Counter-rot. 


Central dip in cr map. 










Dumbbell structure in cr map. 


(g) M9 Ps Db 


Nested IDs (thin and thick) 


Clumps 


Co-rot. 


/z4 > at the ID location. 


(h) M9 Ps Rb 


ID (thin) 


Clumps 


Counter-rot. 




(i)M18Ps Db 


IR+ID (thin) 


Clumps, spiral arms 


Co-rot. 


/z4 > at the ID location. 


G)M18Ps Rb 


ID (thin) 


Irregular clumpy structure 


Counter-rot. 


Central dip in cr map. 


(k)M18Pl Db 


Nested IDs (thin) 


Ring relics 


Co-rot. 




(1)M18P1 Rb 


IR+ID (thin) 


Pseudo-ring 


Counter-rot. 





Columns: (1) Model code. (2) Global structure of the formed IC in each remnant, as derived from the radial and vertical surface density profiles 
shown in Figs. [2ll3l (§ ^3.2. 1113.2.21 ). (3) Substructure found in the face-on surface density maps of Figs.l4ll5l ( ^3.2.31 ). (4) Rotation sense of the IC 
with respect to the global rotation pattern of the remnant (co- or counter-rotating), as obtained in ^3.2.61 (5) Special detectable features present in 
the global ellipticity and position angle (PA) profiles and in the global kinematical maps of the final remnant, pointing to the existence of these 
ICs. All the ICs exhibit abrupt changes of the trend in € and PA profiles at the limits of the ICs, as well as rotation in the center, so none of these 
two features are listed in the Table (see Figs. [T2][T5l in ^3.31 ). 



by the time the satellite nucleus disrupts, but that precession is 
incomplete when the satellite loses its disc. 

The rotation sense of the IC resulting from the disrupted 
satellite material can be deduced from this figure: it rotates in 
the same direction as the main galaxy disc in direct orbits and 
counter-rotates in the retrograde cases. In general, the alignment 
between the formed ICs and the main body of the galaxy is worse 
in retrograde ex periments (as already known, see references in 
lBoumaudll2QQ9l) . Nevertheless, we can conclude that the ICs re- 
sulting in these experiments are highly aligned with the final 
galactic plane in all the cases. 

3.2.6. Rotational support of the ICs 

The ICs formed in our experiments are strongly rotationally- 
supported, with Vmax/cr ~ 2.5 typically. In Fig.[9l we show 
the final velocity and velocity dispersion 2D-maps of models 
M9PsDb and M9PsRb (models g and h in TableO, considering 
all the stars in the remnants, the stars from the primary galaxy, 
those from the satellite disc, and those from the satellite bulge. 
The isophotes of the material considered in each panel are also 
shown, as well as its photometric and kinematic axes. The rota- 
tion velocity maps of both models show that the global velocity 
field of the remnant is governed by the material of the primary 
galaxy. In fact, the existence of a co- (or counter-) rotating IC 
made out of satellite material in the remnant core aflTects neg- 
ligibly to the iso-velocity contours of the global maps (notice 
that the maps considering all the stars and only those from the 
primary galaxy in the remnants are extremely similar in both 
models). Moreover, the existence of these ICs does not aflfect 
appreciably to the velocity dispersion maps either (compare the 
cr maps of each model including and excluding the satellite ma- 
terial). 

Nevertheless, the ICs aflTect noticeably to the orientation of 
the photometric axes with respect to the kinematic ones. As 
shown in the figure, those obtained just considering the primary 
galaxy material are noticeably rotated with respect to the ones 



obtained considering all the stellar content (compare the white 
axes drawn in the two first velocity maps of the direct model). 
This is because direct mino r mergers drive the formation of larg e 
warps in the primary disc (iRoskar et al.ll2010l ISellwoodll2010[) . 
The isophotes in the remnant outskirts of the primary disc ma- 
terial trace these warps, rotating the corresponding photometric 
axes with respect to the orientation (see Fig. [TO]). However, the 
addition of a highly-aligned IC in the centre increases the weight 
of the central regions in the photometry, making the kinematic 
and photometric axes to be more similar (compare the kinematic 
and photometric axes in the first velocity map of the direct model 
in Fig.O. 

On the other hand, a tight alignment between the photomet- 
ric and kinematic axes is observed in retrograde mergers, inde- 
pendently on whether we consider the satellite stars in the maps 
or not. The reason is that retrograde orbits give place to much 
weaker warps than direct ones due to their lower spin-orbit cou- 
pling (compare the velocity maps of the primary disc material in 
both models of the figure). 

The satellite material in the retrograde models counter- 
rotates with respect to the material originally from the primary 
disc (compare the velocity maps of the retrograde model in 
Fig. [9]). However, no counter-rotation is imprinted to the material 
coming from the primary galaxy in any model (but see ^3.2.71) . 
The existence of counter-rotating ICs contribute to rise the ve- 
locity dispersion in the remnant centre, but it does not imprint 
any significant counter-rotation in the global velocity field of the 
remnants. Similar conclusions can be derived from analogous 
maps to the ones shown in Fig.[9]in the other models. 

The findings of § ^3.2.1H3.2.6l shows that minor mergers can 
give place to highly-aligned rotationally- supported ICs out of 
accreted satellite material, without requiring a noticeable re- 
distribution of primary galaxy material through a strong bar or 
relevant dissipative processes (see ®. The rotation sense of 
each IC with respect to the global rotation field of the final rem- 
nant for each model is indicated in Col. 4 of TablelH 
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Fig. 6. Formation of a transitory bar- like distortion in the pri- 
mary disc material during a minor merger. The time evolu- 
tion of the surface density map of the primary disc material of 
model M6PsRs is plotted, using a face-on view of the initial pri- 
mary galaxy (experiment f in Table[3]). Time is shown at the top 
left corner of each panel in simulation units. A rainbow colour 
palette is used to represent different surface density levels in log- 
arithmic scale, with redder colors indicating higher values. All 
physical quantities are given in simulation units. The major axis 
of the oval distortion formed in the remnant disc by the satellite 
impact is marked with a black line in the last panels. It rotates in 
the same direction as the primary disc material (clockwise). The 
orbit and rotation of the accreted satellite is counter-clockwise. 
The life time of this bar-like distortion corresponds to ~ 0.3 Gyr, 
using the scaling indicated in ^ 



3.2.7. Bender diagrams 

Bender diagrams of our remnants have been obtained using the 
procedure introduced by iGonzalez-Garcia et al.l (I2QQ6 ). We ob- 
tain line-of- sight velocity distributions (LOSVD) for the merger 
remnants by choosing a point of view at random for projecting 
the particle distribution. We then derive a surface density map 
and define iso-density contours to fit ellipses deriving values for 
the ellipticity, PA, and the ^4 parameter (the fourth-order Fourier 
coefficient measuring the deviation of the iso-density contours 
from pure ellipses). We place a slit along the major axis of our 
ellipses up to 7? = 2. Given the primary disc scale lengths, our 
mapping reaches well into the region of the disc. We have binned 
the slit in 10 spatial bins and the velocity interval in 50 bins. We 
find the radial projected velocity and the number of particles in 
each bin in velocity for each bin in the slit. In this way we obtain 
a line-of-sight velocity distribution. Finally, we fit the LOSVD 
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Fig. 7. Pseudo-ring structure in the final remnant of model 
M9PsRb, made of stars that initially belonged to the primary disc 
(experiment h in TableO. It coincides with the location of the 
two twin symmetric clumps observed in the stellar satellite ma- 
terial in the final remnant. Both substructures are associated to 
an oval distortion induced by the satellite accretion (see Fig.O. 
The pseudo-ring has been marked with a white dashed circle. 
A rainbow colour palette is used to represent diff'erent surface 
density levels in logarithmic scale, with redder colors indicat- 
ing higher values. All physical quantities are given in simulation 
units. 



by a Gaussian and the residuals by a Gauss-Hermite pol yno- 
mial as given by lvan de r Marel & Franx (1993 ) and Bender e t al.l 
(1994). We repeat this process for each remnant for 90 randomly 
chosen points of view. 

From the fitting procedure we obtain for each LOSVD a 
value for the velocity centroid of the distribution at each spa- 
tial bin along the slit (V), the velocity dispersion (cr), and the 
amplitude of the third Hermite polynomial (h^), which is a mea- 
surement of the skewness of the distribution. 

We show the resulting Bender diagrams for our simulations 
with 1:6 mass-ratio in Fig. [TT] (the results from the other simu- 
lations behave similarly). All models present an anti-correlation 
between the V/cr and parameter. This is something to be ex- 
pected due to the prominence of short-axis tube orbits in the pro- 
genitor discs that are kept relatively undisturbed due to the merg- 
ing event. The large bulge of simulations M6PsDb, M6PsRb, 
M6PlDb, and M6PlRb (labelled from a to d in Table© also 
helps stabilizing these orbits against the eff'ects of the collision- 
less minor merger. These simulations diff'er in the amplitude of 
the skewness parameter or the V/cr as a consequence of the dif- 
ferent orbital parameters of the encounters. 

An interesting diff'erence is to be observed in the two plots 
corresponding to simulations with small bulges in the primary 
disc: M6PsDs and M6PsRs (models e and f in Table[3l). Here 
most of the LOSVD present the anti-correlation usual for discs. 
However for small values of V/cr we find a mild (in the case of 
M6PsDs) or stronger (in M6PsRs) correlation of the two param- 
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Fig. 8. Misalignment between the material initially belonging to 
the satellite bulge and disc with respect to the luminous galac- 
tic plane of the final remnants (diamonds and triangles, respec- 
tively). The letters identifying each model are those used in 
Tables 



and triaxial struct ures in the bulges (see lJungwiert et al] 119971: 
lErwin et"al ] l2QQ3h . 

The existence of weak bar-like or oval distortions in the rem- 
nant center of some models can be deduced from their global 
isophotal profiles. Models M6P1[D/R]b (panels c in Fig.[T2land 
d in Fig. [13]) exhibit a slight maximum in e and a constant PA at 
the radii where the IC presents twin clumps (see the correspond- 
ing panels in Fig.©. Nevertheless, as commented above, none 
of our remnants develop noticeable nuclear bars, a fact that is 
corroborated by the global ellipticity and PA profiles shown in 
Figs.imil 

The models with a small primary bulge exhibit central ellip- 
ticities equal to or greater than that of the outer disc (e ~ 0.3, 
see Fig. [13]), suggesting that the i nner light is dominated by a 
disc or flattened triaxial structure towi n et aP 1200 3 ) . The pro- 
files of these models exhibit ellipticity and PA profiles typical 
of bulge-hdisc galaxies. The ellipticity shows a quick rise in the 
core (bulge) region, a slight decrease in the bulge-disc transi- 
tion region, and a constant value in the disc outer layers, while 
the PA is nearly constant outside the bulge region (see Fig.[T3l). 
Some particular features in the e and PA profile s corroborat ing 
the existence of some of the ICs identified in § ^3.2. 1113. 2. 31 are 
listed in Col. 5 of TableH 



eters. This is due both to the eff'ect of the small bulge that is not 
able to stabilize short-axis tube orbits and to the larger contribu- 
tion of the satellite bulge in the inner parts of these simulations. 
It is interesting to note that the eff'ect of the counter-rotating orbit 
of M6PsRs is signalled by the relatively large correlation signa- 
ture in the central parts of the diagram. 

3.3. Detection of ICs in the global stellar maps of remnants 

The prominence of the primary stars in the final remnant make 
the ICs to be completely masked by their mass distribution, in 
such a way that no direct detection of any IC can be performed 
in either the global surface density maps and profiles of the final 
remnants. It is encouraging though that indirect detection similar 
to that attempted and used in observations proves to be success- 
ful in many cases, as we will show in this section. 

3.3.1 . Ellipticity and position angle profiles of the remnants 

Inner structures in real galaxies imprint noticeable features in 
the ellipticity (e) and PA profiles of the central isophotes of 
the host galaxies (Erwin et al. 20031: iFalcon-Barroso et al ]|200l 
IChemin & Hernandez 2009). In the right panels of Figs. ll2m51 
we have plotted the e and PA profiles of the stellar material in 
all the remnants, assuming an inclined viewing angle (0 = 60°, 
= 20°). The radial extent of the ICs as defined in ^3.2.1! have 
been marked in all the plots (vertical dashed lines). 

The profiles of the remnants with big primary bulges are 
dominated by the light distribution of the pre-existing primary 
bulge in the center, leading to low ellipticities in the remnant 
even at the presence of ICs (e < 0.2, see the cases shown in 
Fig. [121 for example). Nevertheless, the existence of ICs (and 
even their radial extensions) can be deduced just attending to the 
abrupt changes in the trends of the e and PA profiles that appear 
at the transition region between adjacent ICs (for example, be- 
tween two nested IDs or between an ID and the outer remnant 
disc, see Figs. [T2l[T5l) . These trend changes are also observed in 
real galaxies, associated to isophotes twisting produced by bars 



3.3.2. Kinemetric moments of the velocity fields in the 
remnants 

In many real galaxies, the existence of kinematically-decoupled 
components that are not detectable through direct imaging is in- 
ferred from special features in their kinematic maps. In most 
cases, the existence of these ICs is later confirmed by unsharp 
masking or by direct detection of i ts gas component in [QUI] 
or HfS emission-line maps (see, e.g. jFalcon-Barroso et al.ll2QQ3l 
F06). 

Fi gures [T2H 1 5 1 show the 2D-maps of the kinemetric moments 
of the line-of- sight velocity distribution of all the luminous mate- 
rial in each remnant, using edge-on views. The line-of- sight ve- 
locity Vlos, the velocity dispersion cr, and the third- and fourth- 
order coefficients of the Gauss-Hermite expansion hs and /z4 are 
shown for each model. They have been obtained adapting the 
profit routine (originally designed for profile fitting of spectral 
emission-lines by GaussHermite series) to N-body data (iRilfell 
2010). The rotation velocity maps in these figures differ from 
the analogous ones presented in Fig. [9] in the spatial resolution. 
While those presented in Fig. [9] use an uniform spatial binning, 
the ones shown in Figs. [T2l[T5l have higher spatial resolution in 
the center than in the outskirts to improve the S /N of the esti- 
mates of the kinemetric moments in low-density regions (analo- 
gously to what is done in observations, see F06). The central spa- 
tial resolutions in these figures are similar to those achieved by 
current observations in the nearby Universe, adopting the scal- 
ings proposed in ^2] (~ 0.5 kpc). Regions with low S /N have 
been masked in the maps. 

Figures [T2ill5l show that the existence of many ICs could be 
deduced from these maps of kinemetric moments. Many partic- 
ular features present in these maps that corroborate the existence 
of some of the ICs identified in § ^3.2.1113.2.3! are indicated in 
Col. 5 of Table|4j However, we must remark that these features 
are smooth and show large dispersion (specially, at outer radii), 
and hence, they could be detectable only in case of low noise 
levels. Some of these kinematical features include: 

i. Misalignment of kinematic axes. 

A clear misalignment of the inner and outer kinematic axes 
of the galaxy is detected in some cases (see model M6PsRb, 
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panel b in Fig. [121), which can be produced by noticeable disc 
warping (as the case shown in Fig.[TO]). Rotation reaches the 
central regions in all the remnants. The final rotation fields, 
although following a typical spider-like diagram, seem quite 
distorte d in general, as it is characteristic of dissipationless 
models (iJesseit et al.ll2007l) . In general, experiments with lower 
mass ratios or larger pericenters imply a smoother destruction 
of the initial galaxy rotation pattern. 

ii. Stretching of iso-velocity contours in the centre. 

The existence of co-rotating IDs and IRs imprints higher 
rotation in the centre, stretching the iso-velocity contours 
towards the major axis at their location. This makes the angle of 
these contours to be more open at the remnant outskirts than in 
the centre (see, e.g., models M6PsDs and M6PlDb, panels b and 
c in Fig.[T2l). 

iii. Twisting of central iso-velocity contours in retrograde 
models. 

As shown in ^3.2.61 retrograde models give place to 
counter-rotating ICs. In some of them, the contribution of this 
IC to the mass in the remnant core is high enough to imprint a 
strong twisting of the iso-velocity contours by ~ 90° at the core 
region (see model M6PsRb in panel b of Fig. [121 and models 
M18P[s/l]Rb in panels j and 1 of Fig.[l5l). 

iv. S-shaped twists of iso-velocity contours. 

Almost all the cases exhibit S-shaped or integral- sign- 
shaped twisting of the central iso-velocity contours at a certain 
height in the galactic plane. These features are produced by the 
sharp decrease of the amount of stars that contributes to the 
velocity field at a certain spatial position. Therefore, S-shaped 
kinematic twists in our models appear at the spatial locations 
where the mass contribution of the ICs becomes negligible 
(see, e.g., models M6PlRb and MlSPsDb, in panels d and i of 
Figs.[l3land[l4l). 

V. /z4 peaks. 

The majority of the models with big primary bulges show 
positive h4 values in general, peaking at the location of the ICs 
formed at the end of the simulation. In particular, the ring geom- 
etry of the IC resulting in model M6PlRb can be deduced from 
the two well-defined symmetric peaks present in the /z4 maps, 
which trace its perpendicular sections (panel d in Fig.[T3ll. 

vi. Correlations between the kinemetric moments. 

We have found that vlos and correlate at the location 
of the formed IC in our remnants if it exhibits oval or bar-like 
distortions (see models M6Ps[D/R]b and M6PlDb in panels a- 
c of Fig.[l2l and models M9PsRb and MlSPsDb in panels h- 
i of Fig.[T4l). However, when the ICs show no trace of having 
had them, vlos-^3 anti-correlate at the location of the IC, ac- 
cordingly to the discy structure of the majority of the obtained 
ICs (see models MlSPsRb and M18P1[D/R]b in Fig.[l5l). Some 
models also exhibit a mixed behaviour: vlos-^3 correlate at the 
remnant core (at R < 0.3) and anti-correlate at the position of 
the h4 peaks, as occurs in models M6PsDs (panel e in Fig.[T3l) 
and M9PsDb (panel g in Fig. [141). This is probably pointing to an 
origin related to oval distortions and bar-like distortions in these 
components. 

vii. Dumbbell cr structures. 

Dumbbell structures in cr maps point to the existence of 
counter-rotating IDs, which rise cr at their location. Although 
all the retrograde models give place to counter-rotating IDs or 
IRs, dumbbell structures are only observed in models M6PlRb 
and M6PsRs (panels d and f in Fig.[T3b. 

viii. Central cr dips. 
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Fig. 10. Final warped disc of model M9PsDb (experiment g in 
TableO. An edge-on view of the material originally belonging 
to the primary disc in the final remnant of this model is shown. 
A rainbow colour palette is used to represent diff'erent surface 
density levels in logarithmic scale, with redder colors indicat- 
ing higher values. All physical quantities are given in simulation 
units. 



In general, the velocity dispersion is very high in the rem- 
nant centres because the bulge dynamics dominate at those radii. 
However, it tends to decrease at the bulge-to-disc transition 
region due to the contribution of the co-rotating IC. In our 
models, cr dips are detected only in the models of small primary 
bulges (see panels e and f of Fig.[T3l). 

Summarizing, the velocity fields of our remnants exhibit fea- 
tures that reveal the existence of the central ICs formed through 
the satellite accretion, such as the stretching of the iso-velocity 
contours in the centre, S-shaped kinematic twists, and dumbbell 
cr profiles. More inclined views of the remnant smooth out these 
features, making the ICs undetectable in the velocity 2D-maps. 



4. Qualitative comparison to observations 

In this section, we make a qualitative comparison of the prop- 
erties of the ICs obtained in our minor merger simulations with 
real observational cases, to stress the similarities and diff'erences 
among them. 



4.1. Structural comparison 

Fi gure[T6l compares the morphology of some of the ICs formed 
in our experiments to observational e xamples that exhibit sim ilar 
inner structures from the sample of Erwin & Sparke (2003V As 
shown by the figure, the ICs resulting from our minor merger 
experiments are analogous structural and morphologically to 
the ICs hosted by many real spiral galaxies. IDs, nested IDs, 
rings and pseudo-rings, spiral patterns, and even undisrupted 
clumps with irregular spatial distributions as those observed 
in our experiments are found in real spiral galaxies (see also 
Buta & Combes 1996; Erwin & Sparke 2002; Erwin et al. 2 0031 
Erwinii2004: .Perez-Gallego et al.,2010: .SiFchenkoet al...20lTr 
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Star-forming pseudo-rings are substructures u sually associ- 
ated t o gas and starbursts in real galaxies (F06; 'Shapir o et alJ 
1201 (A [Sil'Chenko 2010). However, although our models do 
not consider dissipative components, some remnants exhibit 
stellar ring relics or pseudo-rings in the centre (see models 
M18P1[D/R]b in panels k and 1 of Fig.E]). 

Our remnants have not developed significant nuclear bars 
(although some of them present signs of weak oval or bar-like 
distortions in the center, see ^3.3. lb . Then, the IRs and pseudo- 
rings that appear in some of our models are not directly asso- 
ciated to well-formed bars. Instead, they have formed by the 
coupling of the satellite disruption to the resonances induced 
in the primary disc by the encounter. Cases of IRs and sets 
of concentri c rings not associated to bars are found in real 
galaxies too (jButa & Combeslll996l: iMunoz-Tunon et al" 1[2003; 
iKormendy et al.ll2006h . So. our models point to the resonance- 
related origin of IRs, independently on whether these resonances 
are linked to strong bars or to weak non-axisymmetric perturba- 
tions as the ovals found in some of our models (see ^3.2.41) . 

The stellar IDs formed in our experiments have radial lengths 
ranging from ~ 0.4 to ~ 1 disc scale-length of the final remnant, 
which are equivalent to radial lengths spanning from ~ 360 pc 
to ~ 2.5 kpc, adopting the scaling provided in ^ These values 
are in excellent agreement wit h those obtained for the IDs de- 
tected in rea l galaxies (Erwin & Sparke 20021 ISarzi et al 
IShapiro et al . 2010; Sil'chenko et al. 2011, F06). The IRs and 
pseudo-rings formed in our simulations are embeded in these 
IDs, and thus, exhibit typically lower linear sizes. Adopting the 
scalings of ^ our IRs match pretty well the standard 1 kpc- 
radius IRs that lie betwe en the in ner Lindblad resonances o f 
many disc galaxies (Erwin & Sparke 2002; Shapiro et al.ll2010'). 

Concerning to their vertical distributions, the flared vertical 
structure of most the single IDs and nested IDs resulting in our 
models (see Fig. [3]) is extremely similar to the one derived re- 
cently for the galaxy NCG7217, compose d by two large scale 
nested stellar discs (Sil'chenko et al ]|20ni). 

The existence and extension of many real ICs are usually 
derived from abrupt changes in the trends of the ellipticity 
and PA profiles of the^global galaxy isophotes (lErwin & Sparkel 
[2OO2L 120031: lErwinI l2QQ4l) . As we have shown in ESU aU 
the ICs formed in our models imprint similar features in the 
global pro files of the final remnants, coh erently with observa- 
tions (iPrieto et al.ll200ll: lErwin & Sparke 2003). 

The low misalignment of the ICs obtained in the simulations 
with respect to the global galaxy plane (see ^3. 2. 5b is in excel- 
lent agreement with observations: F06 report misalignments of 
typically < 20° for th e ICs harboured by Sa-Sb galaxies, and 
'Krainovic et al.' ('2011) find that 90% of the ICs in a sample of 
260 early-type galaxies can be considered aligned to better than 
5 deg. Nevertheless, we can not discard the possibility that this 
may be facilitated by the moderate inclinations of our simulated 
orbits. However, the quick decay of the satellite orbit to the pri- 
mary disc plane observed in all the models suggests that more 
inclined orbital configurations could result into aligned ICs too, 
and that extreme orbital inclinations would be required to pro- 
duce noticeably misaligned ICs. 

Our simulations also demonstrate the feasibility of the sce- 
nario proposed by F06 to explain the low misalignment of the 
ICs observed in Sa-Sb galaxies as compared to those in E-SO's 
(usually above 60°, see lSarzi et al ]|2006). According to these au- 
thors, if the ICs have resulted from minor mergers, the morphol- 
ogy of the primary galaxy in the encounter should play a crucial 
role in determining its alignment. In our simulations, the evo- 
lution of the satellite orbit to the primary disc plane prior to its 



disruption is produced by the low spheroidality of the primary 
galaxy potential, which establishes a privileged plane a priori 
(the primary disc plane). Spheroidal potentials (as those of E- 
SO galaxies) do not have any privileged direction, and thus they 
do not produce this eff ect in the orbits of the accreted satellites 
(iDi Matteo et al.|[2008h . Therefore, our simulations suggest that, 
if the majority of ICs derive from minor mergers, we should ex- 
pect to detect more aligned ICs in disc galaxies (as Sa-Sb's) than 
in spheroidal ones (E-SO's), coherently with observations and 
with F06 arguments. Moreover, F06 propose that the existence 
of large gas amounts involved in the merger must contribute to 
this alignment. Notice that, although the previous sentence must 
be true, our models prove that gas is not strictly necessary to 
obtain high co-planar ICs in minor mergers. 

4.2. Kinematical comparison 

ICs disturb the velocity field of their host galaxies, giving place 
to particular kinematical features in the global maps of the 
galaxy pointing to their existence. Our remnants exhibit many 
of these features, such as disturbed iso- velocity contours in the 
center, noticeable S- shaped kinematic twists at the limiting radii 
of the ICs, cr peaks, dumbbell-like cr structures associated to 
counter-rotating IDs, and stretching of iso-velocity contours at 
the presence of co-rotating ICs (see ^3.2.61) , similarly to w hat it 
is observed in real gal axies harbouring ICs (see F06; Sarzi et al.l 
I2OO6I: iKrainovic et aD l2008). 

The counter-rotating ICs formed in our models are so embe- 
ded in the primary disc stellar material that, in general, they do 
not im print n oticeable counter-rotation at the centre of the rem- 
nant ( ^3.3.21) . The masses of the ICs are too low to affect signif- 
icantly to the central dynamics of the remnants, which are ba- 
sically dominated by the host rotating disc instead (see ^3.2.61) , 
although some counter-rotation becomes detectable at specific 
viewing angles (see model M6PsRs in panel f of Fig.fTTl). 
Therefore, the presence of large disc components in galaxies 
makes the observational detection of counter-rotating ICs dif- 
ficult. This can explain why the observational samples of early- 
type spirals exhibit an apparent absence of counter-rotating fea- 
tures as compa red to the samples of E-SO's, w hich lack of rele- 
vant discs (see iKannappan & Fabricantll200ll F06). In fact, the 
frequency of counter-rotating ICs in E-SO's and in Sa-Sb's could 
be similar, the discs of Sa-Sb's being responsible of masking 
them. 

Our collisionless models reproduce several trends and cor- 
relations between the kinemetric moments of real galaxies with 
ICs, although they present differences as well. In general, our 
models show /z4 > and h4 peaks at the location of the ICS, 
in excellent a greement with observ ations of ICs in early-type 
galaxies (^F06: lKrajnovic et aDl2008l) . However, the Vlos-^3 cor- 
relation found in our simulations reproduce the observational 
trends partially. Observations report Vlos-^3 anti-correlation at 
the location of co-rotating ICs in Sa-Sb's, and correlation if 
the ICs counter-rotate (F06). Nevertheless, our models tend to 
exhibit Vlos-^3 anti-correlation, independently on the rotation 
sense of the IC with respect to the final galaxy (see Figs. [12]- 
[T5l) . Although some retrograde models can present correlation 
for particular viewing angles (see model M6PsRs, panel f in 
Fig.[TT]), in general we find anti-correlation (see the rest of 
Bender diagrams in Fig.fTTl). 

Our models do not include the effects of gas dynamics and 
star formation, which obviously affect to the kinematical struc- 
ture of the final remnant (see Collisionless remnants are 
known to trigger box orbits, which move Vlos and hs to cor- 
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relate (see, e.g. jNaab & Burkertll200lh . However, the inclusion 
of gas in the simulation has the opposite effect: it suppresses 
ther n and make Vlqs and to anti-correlate at the remnant cen- 
tre (iBendo & Barnesll200a iNaab & Burkertll2Q0ll iJesseit et al.l 
[2QQ7ir Then, our collisionless simulations behave atypically in 
this sense, as they tend to exhibi t Vlqs -^3 anti-correlati o n. Thi s 
trend was already obtained by "Gonzalez- Garcia et al.l (l2QQ6l) . 
who showed that such anti-correlations can be kept in a colli- 
sionless merger simulation whenever a central bulge allows the 
discs to retain some of their original angular momentum during 
the merger, making short-axis tube orbits to be still present in 
the final remnant. 

In our models, Vlos-^3 correlation seems to be more re- 
lated to the existence of bar-like or oval distortions than to 
counter-rotation ( S3. 3. 2b . The behaviour observed in our rem- 
nants can be reconciled with the observational trends if we con- 
sider that retrograde mergers give place to longer-lasting bar- 
like instabilities (although weaker) than their direct analogues 
(iBames & Hemquisil 19961) . Therefore, we should expect to find 
statistically more bar-like distortions and ovals in the ICs result- 
ing from retrograde mergers than in those coming from a direct 
encounter. Assuming that Vlos-^3 correlation is related to the 
existence of non-axisymmetric distortions in the galaxy disc, the 
previous result establishes a correspondence between Vlos-^s 
correlations and counter-rotating ICs, coherently with F06 re- 
sults. This scenario is supported by the diff'erent Vlos-^b trends 
exhibited by Sa-Sb and E-SO galaxies. While ICs hosted by 
the later ones are alw ays associated to Vlos-^3 anti-correlations 
(iBender et al.lll994l) . the ICs in early-type discs exhibit corre- 
lations only if the IC is counter-rotating (F06). This can be ex- 
plained considering that E-SO's are more efficient inhibiting bars 
and ovals than Sa-Sb's (because they have larger bulges and neg- 
ligible disc components), a fact that should make the ICs in these 
galaxies to exhibit Vlos-^s anti-correlations in general (coher- 
ently with observations). 

In conclusion, minor mergers could account for the existence 
of many stellar dynamically-cold ICs in spiral galaxies and even 
in E-SO's. The present models prove that they can give place 
to ICs with geometrical, structural, and kinematical properties 
similar to those observed in real galaxies. 

5. Model limitations 

We have analysed the role of minor mergers in the formation 
of dynamically-cold thin stellar ICs, without accounting for gas 
and star formation eff'ects. Obviously, the inclusion of dissipative 
components in the models would not just provide the formed IC 
with an additional recent stellar population or gas component, 
but it could aff'ect noticeably its final structure and kinematics. 

The formation of IDs and IRs are associated to the redistri- 
bution of angular momentum in the remnant disc. In this sense, 
gas components are expected to contribute noticeably to this 
re-distribution, mainly during disc distortions. Nevertheless, re- 
cent studies have demonstrated that dissipative components (al- 
though relevant) are not essential or decisive in the formation of 
kinematically-decoupled ICs through major mergers. Gas causes 
the remnants to appear more round and axisymmetric, wiping 
out small kinematical misalignments more easily, but the result- 
ing IC conserves mos t of its structural p roperties as compared to 
collisionless models ( IJesseit et al.ll2QQ7l) . 

Besides, although star formation triggers the formation of 
young stars in the IC in simulations of major mergers, the re- 
sulting IC is still com posed by a relevant old stellar component 
(iDi Matteo et al.ll2008i) . If the formation of an ID or IR depended 



basically on merger-induced gas inflows to the galaxy cen- 
tre and on the subsequent star formation (Barne s & HemquistI 
Il996; Bournaud et al. 2005b), these substructures should be 
bluer or, at least, younger than the surrounding bulge compo- 
nent. However, this is not the case, as IDs exhi bit very simi- 
lar colors to those of their host bulges usually (see iMorelli et al.l 
120041: iPeletie" et al.l l2007l) . Therefore, although gas and star for- 
mation eff'ects must have been relevant in the formation of IDs 
and IRs, they may not be essential for it in many cases. This 
is also supported by the significant old underlying stellar com- 
ponent detected in most IDs and IR s (Buta & Purcefl 1998; 
Buta et al.' '1998; van den Bosch et al.l [T998I: iKrainovic & Jaff'j 
2004; MorelH et al. 200^^ 

The flattened structure of IDs and IRs has been tradition- 
ally interpreted as a sign of the essenti al role played by gas in 
their formation (ICappellari et al 1 120071) . but our models prove 
that central thin rotationally-supportedlCs can result from satel- 
lite disruption without supplying gas to the remnant centre. 
However, this is a simplified picture of the reality, as the major- 
ity of the observed IDs and IRs contain recent (or on-going) star 
formation, dust, and gas, a fact that clearly points to the tight re- 
lation between dissipative processes and their buildup (see, e.g. , 
Barnes & H emaui st* '1996^ 'va n den Bosch & EmsellemI 119981: 
MorelH et al. 2004; Kormend vetaD l2005l: ICappeflari elaP 
2007: iPeletier et al.i»200 7. F06). In fact, we could expect that the 
inclusion of gas and star formation in the models would made 
the formed ICs more detectable, accounting for the star forma- 
tion t hat minor mergers u sually trigger in the centre of galaxy 
discs (iKaviraj et al.ll2009b . So, we intend to re-run these simu- 
lations including gas and star formation processes in the near 
future, to determine the eff'ects of dissipative processes in the 
formation of ICs. 



6. Discussion 

The present models provide a novel insight into the buildup of 
dynamically-cold ICs through minor mergers, because the re- 
sulting IDs and IRs do not derive from primary disc material 
or from a gaseous component as in other studies, but from the 
disrupted stellar satellite material. Traditionally, minor merg- 
ers have been considered as secondary drivers in the forma- 
tion of IDs or IRs in galaxies, in the sense that the rotationally- 
supported ICs were thought to come from the bars triggered in 
the discs by the encounters, not from the minor merger them- 
selves. The accreted satellite material did not ended in these ICs 
except if it reached the centre undisrupted (see ^T])- Therefore, 
minor mergers have been considered just as the agents inducing 
bars, but the bars were the processes responsible for the buildup 
of the ICs. 

Our models show that both processes (the induced disc res- 
onances and the accretion of material external to the galaxy) are 
extremely connected, as the resulting IC is made out of satellite 
material and exhibits an origin related to transitory non axisym- 
metric distortions of the galaxy disc at the same time (§ S3.U 
13.21) . The diff'erence between our models and previous ones from 
the literature lies basically in the satellite characteristics, as we 
have found that very diff'erent orbital configurations give place to 
similar ICs (see Table|4]). In our models, the discy structure and 
the realistic density contrast of the satellites with respect to the 
primary galaxies make them to be sensitive to disruption, but re- 
sistant enough to reach the remnant centr e (in cont rast with pre- 
vious models, see Aguerri et al. 2001; M oster et al . 2010). This 
implies that, if a low-density satellite (as a dS) were accreted 
by a galaxy, the existence of a prominent bulge in the primary 
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galaxy would induce the complete disruption of the satellite and 
the formation of an IC out of it. However, if the accreted satel- 
lite were dense (such as a dSph or dE) or if the primary galaxy 
had a small bulge, the undisrupted satellite core would deposit 
in the centre without disrupting, increasing the spatial density in 
the galaxy core. 

The bar-related origin of numerous bulges and dynami cally- 
cold ICs in gala xies is i ndispu t able ( Martinez- Valpuesta et al] 
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20061: Berentze n et al.1 I2Q07I: iRoniero-Gomez et al.l 
Athanassoula et al.1 l2009allbr 120101: iBaglev et al.l 
Fisher & Drorv 201 Oh . but it is also true that more than 
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one third of them in Sa-Sb galaxies are apparently unrelated to 
bars (see references in Considering the relevanc e that minor 
merg e rs seem to have had in galax y evolution (see I Jo gee et aP 
120091: iLopez -Saniua n et al.1 l2010all b'). it is probable that even 
many bars may have been induced by minor mergers in turn. 
So, the present models provide a feasible explanation to the 
existence of old, pure stellar IDs and IRs i n unbarred galaxies 
(specially to the co unter-rotating cases, see iJesseit et al .1120071: 
iMorelH et al.ll2010h . 

The mixing of material observed in the models (due to the 
merger and/or to the triggered oval distortions) indicates that 
it would be very difficult to prove the minor merger origin of 
the resulting ICs by disentangling the stellar population with an 
external origin from the underlying stars. However, the minor 
merger origin of the ICs in some galaxies is obvious, mainly 
if the IC is counter-rotating or it is harboured by a non-relaxed 
host ( see'Hay nes et al.ll2000l:lGutierrez et al.ll2002l: IShapiro et all 
[20TQl: [Sirchenk o et al.ll201lh . 

The key role of minor mergers in the growth of bulges 
of spiral galaxies is evident, as relics of disrupted satellites 
and on-going minor mergers are frequently found in nearby 
galaxies, indicating that these p rocesses are extremel y common 
(iMartmez-Delgado et al.l 120071 I2008L I2009L 120101: iKnieni^ 



2010). Therefore, the present models suggest that the majority of 
the ICs found in spiral galaxies could have had a minor-merger 
related origin (independently on whether the minor merger trig- 
gers a noticeable bar or not), and that the role of minor mergers 
in the formation of ICs may have been much more complex than 
just bar triggering, as traditionally assumed. 

7. Summary and conclusions 

We have investigated the capability of minor mergers to trigger 
the formation of IDs and IRs in spiral galaxies through collision- 
less N-body simulations. We have extended the simulations of 
minor mergers onto disc galaxies presented in EM06, sampling 
a wider parameter space of initial conditions. Different orbits 
and mass ratios have been considered, as well as two different 
models for the primary disc galaxy (Sab or Sc galaxy). 

All the simulated minor mergers have developed thin 
rotationally- supported ICs out of disrupted stellar satellite ma- 
terial, with scale lengths analogous to those observed in real IDs 
and IRs. The resulting ICs are highly aligned to the main galac- 
tic plane of the remnant, as the original non- spheroidal potential 
of the primary galaxy makes the satellite orbit to evolve to its 
privileged plane prior to disruption. This fact provides a possi- 
ble explanation for the low misalignment observed in the ICs 
found in Sa-Sb galaxies as compared to those in E-SO's. No rel- 
evant counterparts of these dynamically-cold ICs are obtained in 
the remnant material originally belonging to the primary galaxy. 

The geometrical analysis of these ICs reveals a wide mor- 
phological zoo of ICs, similar to the observed ones in Sa-Sc 
galaxies through direct imaging or through unsharp masking. 



such as IDs, IRs, pseudo-rings, nested IDs, spiral structure, and 
different combinations of them. No noticeable bars have been 
formed either in the primary disc or in the IC made out of ac- 
creted satellite material. The structural and kinematic properties 
of these ICs are analogous to those observed in real galaxies as 
well. The existence of these ICs cannot be derived directly from 
global surface density maps and profiles, but it can be deduced 
from the characteristic features that they imprint to the isophotal 
profiles and kinematic maps of the final remnants, analogously 
to many observational cases. 

Key points to explain the formation of these IDs and IRs in 
our simulations without requiring significant bars or dissipative 
components are the structure and density of the satellites, as well 
as the existence of a prominent bulge in the primary galaxy. The 
realistic satellite-to-primary galaxy density ratios make the satel- 
lites to be more sensible to orbital circularization and disruption 
than those used in previous simulations. 

Combined with the disc resonances induced by the en- 
counter, these three processes (satellite disruption, orbital cir- 
cularization, and coupling of disruption with merger- triggered 
resonances in the disc) give place to highly aligned rotationally- 
supported ICs in the remnants centre. The existence of big bulges 
in the primary galaxies and long-lasting decaying orbits ensure 
a more efficient satellite disruption, producing thinner ICs. This 
implies that, if a low-density satellite (such as a dS) were ac- 
creted by a galaxy with a prominent bulge, it would result in 
the complete disruption of the satellite and the formation of a 
rotationally- supported flat IC. However, if the galaxy accreted 
were a high-density satellite (such as a dSph or dE) or if the pri- 
mary galaxy had a small bulge, the undisrupted satellite core is 
expected to sink to the galaxy centre without disrupting com- 
pletely, contributing to the formation of a central bulge. In this 
sense, minor mergers could account for the existence of old, 
pure stellar IDs and IRs in many unbarred galaxies (specially 
the counter-rotating ones). 

Traditionally, minor mergers have been considered just as 
secondary agents in the formation of dynamically-cold ICs in 
galaxies, only responsible for inducing the bars in the galaxy 
discs that finally give place to these ICs. Our models suggest 
that the majority of the ICs found in spiral galaxies could have 
had a minor-merger related origin (independently on whether the 
minor merger triggers a bar or not), and that the role of minor 
mergers in the formation of ICs may have been much more com- 
plex than just bar triggering. In conclusion, the present models 
prove that minor mergers are an extremely efficient mechanism 
to form rotationally- supported stellar ICs in spiral galaxies, nei- 
ther requiring strong dissipation nor the development of strong 
bars. 
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Fig. 2. Radial surface brightness profiles of disrupted satellite material in the final remnants of all the models. The plotted physical 
magnitudes are provided in simulation units. The letters identifying each panel are those used in Tabled Circles: Considering all 
the stars originally belonging to the satellite. Red dashed lines: Satellite bulge stars. Blue dotted lines: Satellite disc stars. Black 
solid straight lines: Exponential fits to the radial density profiles considering all the satellite stars. The data points plotted in green 
have been excluded from the fits. In the panels corresponding to remnants that are better fitted with two exponential discs, the two 
scale-lengths are shown in the frame. When just one exponential disc provides a better fit to the global satellite stellar profile than 
two nested ones, only one scale-length is provided. The root mean square of the global fit is also shown in each remnant. Orange 
light solid curved lines: Total fit to the radial surface density profiles in those cases that are better fitted by two exponential discs. 
Vertical dashed lines: Extent of each IC characterized by each fitted exponential radial profile, as defined in the text. The radial 
extension of the ICs is remarked using a shaded background with diff'erent color for the diff'erent ICs identified in each model (green 
for the innermost ones and orange for the outer ones). Sub-panels: Relative residuals of the fits as a function of the radial position. 
The models with small primary bulges obey a diff'erent scaling to the rest of models (see ^2]). Therefore, the simulation units in these 
two experiments diff'er from those in the others, as observed in Figs.|2]-|4] 
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a) M6PsDb: Nested IDs b) M6PsRb: Nested IDs 
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c) M6PIDb: Nested IDs d) MGPIRb: IR embeded In ID 
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e) M6PsDs: Bulge + ID f) M6PsRs: Bulge + ID 
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Fig. 4. Morphology of the ICs made out of disrupted satellite stellar material resulting in models with mass ratio 1:6 (models a-f in 
TableO. Surface density maps of the stars originally from the satellite in the final remnants are presented, using a face-on and an 
edge-on view (left and right columns in each model, respectively). A rainbow colour palette is used to represent diff'erent surface 
density levels in logarithmic scale, with redder colors indicating higher values. The levels of the color scale diff'er from panel to 
panel, as they have been set automatically to ensure an adequate sampling of the dynamical range of the surface density shown in 
each map. Spatial scales in both axes are pro vided in simulation units. The radial extent of each independent IC (i.e., characterized 
by an unique exponential profile, see S3. 2. lb are marked with dashed lines in each map. 
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g) M9PsDb: Nested IDs 



h) M9PsRb: ID 
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i) MlSPsDb: IR embeded in ID 



j) MlSPsRb: ID 
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k) MlSPIDb: Nested IDs 



I) MlSPIRb: IR embeded in ID 
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Fig. 5. Morphology of the ICs resulting in models with mass ratio 1:9 and 1:18 out of disrupted satellite stellar material (models g-1 
in TableO. See caption of Fig.|4l 
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g) Model M9PsDb 



h) Model M9PsRb 
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Fig. 9. Line-of-sight rotation velocity, velocity dispersion, and signal-to-noise ratio (S /N) maps of the stellar material coming from different original components in the final 
remnants of models M9PsDb and M9PsRb, using an edge-on view of the remnants (experiments g and h in TableO. The same levels have been used for the color palette of the 
velocity maps of the diff'erent components in each model and in the velocity dispersion maps to remark the high rotational support of the ICs resulting from satellite material 
(consult the bars on the top of the corresponding columns). Levels are distributed linearly. The maximum level in the bars corresponding to the S /N maps show four diff'erent 
values, corresponding to each one of the frames below, respectively. All physical magnitudes are provided in simulation units. First row: Considering all the stars in the remnant. 
Second row: For the stars initially belonging to the primary galaxy. Third row: For the stars originally in the satellite disc. Fourth row: For the stars originally in the satellite 
bulge. Contours: Surface iso-density contours in the final remnant of the material considered in each map, just as reference. White straight lines: Photometric axes of the material 
considered in each panel. Black straight lines: Kinematical axes of the material considered in each panel. 



Eliche-Moral et al.: A minor merger origin for stellar inner discs and rings in spiral galaxies 




-3 -1 '% 12 3 -3-2-10123 -3 -2 -1 D 1 2 3 

V/ff W 

Fig. 11. Bender diagrams of the remnants resulting from the minor merger experiments run with a mass ratio 1:6. The amplitude of 
the third Hermite polynomial h^, of the final remnant is plotted against its Vjcr for 90 randomly chosen points of view. 
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a) MfiPsDb: Nested IDs 
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b) MGPsRb: Nested IDs 




c> M6P1Db: Nested IDs 
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Fig. 12. Photometrical and kinematical features imprinted by the formed ICs in the global stellar maps of the final remnants of 
models M6Ps[D/R]b and M6PlDb (models a to c in TableO. Left panels: Ellipticity and PA isophotal profiles of all the stars in the 
remnants, using an inclined view (6 = 60°, = 20°). The extent of the ICs as defined in ^3.2. H is marked with vertical dashed lines 
in each case for reference (accordingly to Fig.|2). Right panels: 2D-maps of the kinemetric moments of the line-of-sight velocity 
distribution of all the stars in the remnant using an edge-on view (Vlos, cr, h^, h/C). The levels are distributed linearly (simulation 
units). The photometric axes defined by the material in the core region of the galaxy are also plotted {straight lines). The isophotes 
of the stellar material in the remnant originally belonging to the satellite are overplotted just as reference {contours). 
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d) MGPIRb: IR + ID 
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Fig. 13. Photometrical and kinematical features imprinted by the formed ICs in the global stellar maps of the final remnants of 
models M6PlRb and M6Ps[D/R]s (models d to f in Table©. See caption of Fig.IHl 
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g) M9PsDb: Nested IDs v,. 
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i) MlSPsDb: IR + ID 
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Fig. 14. Photometrical and kinematical features imprinted by the formed ICs in the global stellar maps of the final remnants of 
models M9Ps[D/R]b and MlSPsDb (models g to i in Table©. See caption of Fig.fHl 
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j) MlSPsRb: ID 
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Fig. 15. Photometrical and kinematical features imprinted by the formed ICs in the global stellar maps of the final remnants of 
models MlSPsRb and M18P1[D/R]b (models j to 1 in Table[3]). See caption of Fig.[T2l 
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MODELLED ICs REAL ICs MODELLED ICs REAL ICs 

b) M6PsRb NGC 4665 d) M6PIRb NGC 2273 




g) M9PsDb NGC 5377 h) M9PsRb NGC 6654 




i) MlSPsDb NGC 4314 j) MlSPsRb NGC 4691 




k) MlSPIDb NGC 4245 I) MlSPIRb NGC 718 




Fig. 16. Comparison of some ICs obtained in our minor merger experiments (first and third columns in the figure) to real ob- 
servational examples with similar morphologies (second and fourth columns, respectively). The surface density maps of the ICs 
generated in our models are taken from Figs. |4]|5] and are in simula tion units (consu l t cap tions there). The observational examples 
are taken from the sample of spiral galaxies with ICs developed by Erwin & Sparke (2003 ). Coloured maps of real galaxies corre- 
spond to B - Rot V - I color maps of these galaxies, while their grey-scale maps represent unsharp masks in V, Ror H bands (see 
lErwin & Sparkel (l2003l) for a detailed description of each frame). The global inclination of galaxy in the observational cases and the 
spatial scale of each postage stamp are indicated in each frame. In images of real galaxies. North is up and East is left, except in 
NGC 4665, which has been rotated 90° clockwise to emphasize the similarity with the IC resulting in experiment M6PsRb (model 
b). For a comparison o f the scaled sizes of the ICs resulting in our models with those exhibited by real ICs, consult ^4] Material 
m)m lErwin & Sparkel in this figure is reproduced by permission of the American Astronomical Society (AAS) and of the 

original authors. 



